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ABSTRACT 

With an aim of investigating the nature of evolution-induced mixing in the enve- 
i— i lope of evolved intermediate-mass stars, we carried out an extensive spectroscopic 

study for 12 Cepheid variables of various pulsation periods (~ 2-16 days) to deter- 
j^yy mine the photospheric abundances of C, N, O, and Na, which are the key elements 

for investigating how the H-burning products are salvaged from the interior, based 
on 122 high-dispersion echelle spectra {y 10 per target) of wide wavelength coverage 
collected at Bohyunsan Astronomical Observatory. Having established the relevant at- 
mospheric parameters corresponding to each phase spectroscopically from the equiv- 
alent widths of Fe I and Fe II lines, we derived C, N, O, and Na abundances from C I 
7111/7113/7115/7116/7119, O I 6155-8, N I 8680/8683/8686, and Na I 6154/6161 lines 
by using the spectrum- synthesis fitting technique, while taking into account the non- 
LTE effect. The resulting abundances of these elements for 12 program stars turned out 
to show remarkably small star-to-star dispersions (< 0.1-0.2dex) without any signifi- 
^ cant dependence upon the pulsation period: near-solar Fe ([Fe/H] ~ 0.0), moderately 

CO underabundant C ([C/H] ~ -0.3), appreciably overabundant N ([N/H] ~ +0.4-0.5), 

ON and mildly supersolar Na ([Na/H] ~ +0.2). We conclude the following implications 

from these observational facts: (1) These CNO abundance trends can be interpreted 
mainly as due to the canonical dredge-up of CN-cycled material, while any significant 
£C) non-canonical deep mixing of ON-cycled gas is ruled out (though only a slight mixing 

may still be possible). (2) The mild but definite overabundance of Na suggests that 
CO the NeNa-cycle product is also dredged up. (3) The extent of mixing-induced pecu- 

liarities in the envelope of Cepheid variables is essentially independent on the absolute 
magnitude; i.e., also on the stellar mass. 

•i-H 
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1 INTRODUCTION 

As a star is evolved off the main sequence after exhaustion of 
hydrogen fuels in the core, it increases its radius while the 
surface temperature drops down, and the deep convection 
zone is developed, by which part of the nuclear-processed 
products in the core may be salvaged and mixed into the 
outer envelope and abundance peculiarities may be observed 

* Based on observations carried out at Bohyunsan Astronomical 
Observatory of Korean Astronomy and Space Science Institute 
(KASI). 

f Large data of electronic tables are provided as supplementary 
materials. 
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for some specific elements. By making use of this fact, we can 
study the physical process in the invisible interior of stars by 
comparing spectroscopically determined surface abundances 
with theoretical expectations from stellar evolution calcula- 
tions. 

While considerable progress has been made so far in this 
field and important observational characteristics are known 
to be successfully explained by theoretical calculations, dis- 
crepancies between the empirically established surface abun- 
dances and the prediction from the standard theory are oc- 
casionally seen. One of such problems concerns the oxygen 
abundances in F-G supergiants, which is namely the dis- 
agreement between the theoretical prediction (almost nor- 
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mal O) and an apparent deficit often suggested from obser- 
vations. 

In the envelope of such intermediate-mass supergiants, 
the H-burning (CNO cycle) products are mixed to alter 
the surface abundances. According to the canonical stellar 
evolution calculations (e.g., Lejeune & Schaerer 2001), it 
is essentially the CN-cycled (C— reaction) material that 
is dredged up, while the product of ON-cycle (O— »N reac- 
tion; occurring in deeper region of higher T) is unlikely to 
cause any significant abundance change because mixing is 
not expected to substantially penetrate into such a deep 
layer; thus the predicted surface abundances are character- 
ized by a deficit in C as well as an enhancement in N, while 
O is practically unaffected (if only slightly decreased). 

On the observational side, however, Luck & Lambert 
(1981, 1985) reported based on the CNO abundance results 
of FGK supergiants (including Cepheid variables) that oxy- 
gen is mildly deficient relative to the Sun by ~ 0.2-0.4 dex, 
while C is underabundant by ~ 0.2-1.0 dex and nitrogen 
is overabundant by ~ 0.2-0.8 dex (both with a large star- 
to-star dispersion). In such a deficiency of O is real, a non- 
canonical deep mixing may exist that dredges the ON-cycle 
product up into the envelope. Yet, Luck & Lambert (1985) 
did not stick to this solution and discussed several possibil- 
ities (e.g., errors in abundance determinations, peculiarities 
in solar abundances, etc.). Their tentative conclusion was 
that their absolute abundance determinations may not nec- 
essarily be trustworthy, pointing out that the trends of more 
reliable C/O ratios (from forbidden lines) as well as C/N ra- 
tios (from permitted lines) could be explained by canonical 
dredge-up of CN-cycled material (though its amount may be 
more severe than expected). Thus, the problem remained as 
an open question. 

Meanwhile, Takeda and Takada-Hidai published a series 
of papers in 1990s (Takeda & Takada-Hidai 1994, 1995, 1998, 
2000), where the abundance peculiarities of supergiants in 
terms of Na, N, O, and C were investigated by taking into 
account the non-LTE effect, though they did not pay much 
attention to determination of atmospheric parameters (ap- 
propriate values corresponding to the spectral type and the 
luminosity class were simply assumed). One of the important 
results they first learned was the confirmation of Na enrich- 
ment (Takeda & Takada-Hidai 1994), which is regarded as 
an evidence for the dredge-up of NeNa-cycle product (cf. 
Sasselov 1986, Lambert 1992). Since such a peculiarity in 
Na was not expected from canonical stellar evolution cal- 
culations, it implied an instructive fact that the standard 
theory was not necessarily sufficient. 

Regarding CNO, Takeda & Takada-Hidai (1998) con- 
firmed the mild underabundance of O by ~ 0.3 dex, interest- 
ingly just like Luck & Lambert (1985) reported, along with 
the enrichment of N (Takeda & Takada-Hidai 1995) as well 
as the deficit of C (Takeda & Takadai-Hidai 2000; though 
only for late-B through late- A supergiants). But Takeda & 
Takada-Hidai (1998) did not associate this result of subso- 
lar O with non-canonical mixing of ON-cycled material but 
attributed it to the general tendency of apparently subso- 
lar CNO observed in young stars such as B-type stars (e.g., 
Nissen 1993), which was generally accepted (though puz- 
zled) at that time. However, an extensive non-LTE study on 
the oxygen abundances of B-type stars in comparison with 
the Sun recently carried out by Takeda et al. (2010) lead 



to the conclusion that B stars have almost the solar com- 
position in terms of O without any significant difference. 
If so, these previous studies may point to the conclusion 
that O is apparently underabundant in supergiants while 
primordial O must have been almost normal when they were 
born. Does this imply that oxygen abundance in the enve- 
lope of intermediate-mass stars really got decreased during 
the course of evolution by a non-standard deep mixing of 
ON-cycled product? 

In the meantime, Kovtyukh & Andrievsky (1999) sug- 
gested a possible key to the solution of this oxygen problem. 
According to them, the atmospheric parameters adopted by 
previous investigations may not have been appropriate, since 
spectroscopically determined parameters based on Fe I and 
Fe II lines tend to suffer to appreciable errors because of 
the non-LTE overionization affect (particularly important 
for stronger Fe I lines). As a method to alleviate this prob- 
lem, they proposed to principally invoke Fe n lines (consid- 
ered to be hardly affected by any non-LTE effect) of vari- 
ous strengths for determination of spectroscopic parameters 
while using only "fairly weak" Fe I lines to define the refer- 
ence abundance derived from neutral Fe lines (which is to 
be made consistent with that from once-ionized Fe lines for 
the requirement of ionization equilibrium) as the weak-line 
limit under the assumption that very weak deep- forming Fe I 
lines are formed nearly in LTE. They tested this method 
by applying it to analyzing the spectra of 5 Cep at seven 
different phases, and obtained quite consistent results be- 
tween the such derived spectroscopic gravity and the physi- 
cal gravity (determined from mass and radius), which even- 
tually solved the long-standing gravity-discrepancy problem 
in Cepheids (e.g., Luck & Lambert 1985). Then, as an im- 
portant consequence, the resulting oxygen abundance with 
the new parameters (higher logg) turned out to be nearly 
normal, in contrast to the case of old parameters (lower log g; 
derived by the conventional approach) where an appreciable 
underabundance of oxygen was obtained. So, if their claim 
is correct, the oxygen problem in F-G supergiants may be 
nothing but an apparent effect due to an inadequate choice 
of atmospheric parameters; and the "actual" surface CNO 
abundances would be characterized by low C, high N, and 
nearly normal O, which may be reasonably explained by the 
theory of canonical mixing (i.e., essentially the dredge- up of 
CN-cycle product). 

Motivated by the argument of Kovtyukh & Andrievsky 
(1999), we decided to carry out an extensive spectroscopic 
analysis on the abundance peculiarities of C, N, O, and 
Na (the key elements for studying the mixing of H-burning 
products) for Cepheid variables of various pulsation periods 
based on high-dispersion spectra taken at different phases, 
in order to see whether their conclusion can be generally 
confirmed for evolved intermediate-mass stars of diversified 
properties, which is the purpose of this study. The reason 
why we have chosen Cepheids (instead of non-variable su- 
pergiants) is because (1) we can directly estimate the un- 
certainties in abundance determinations by comparing the 
results from spectra of different phases, and (2) the stellar 
parameters (such as absolute magnitude, mass, radius, etc.) 
are well known in advance by making use of the period- 
luminosity relation. 

What we want to check in this investigation are the fol- 
lowing points : 
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— Are we able to derive reliable atmospheric parameters 
spectroscopically based on Fe I and Fe n lines by the way 
proposed by Kovtyukh & Andrievsky (1999)? Can we ob- 
serve an agreement between the spectroscopic gravity and 
the dynamical gravity? 

— Are the abundances derived from spectra corresponding 
to various phases reasonably consistent with each other? 

— What about the tendency of C, N, O, and Na abun- 
dances in view of the specific trends resulting from mixing 
of H-burning products? In particular, how is the abundance 
of oxygen? Can we confirm that the sum of C+N+O is con- 
served? 

— Is there any meaningful dependence in the extent of abun- 
dance peculiarities upon the stellar mass (or luminosity di- 
rectly related to the pulsation period), such as suggested 
previously (e.g., Takeda & Takada-Hidai 1994) for the case 
of Na? 

It may be worth noting here that extensive spectro- 
scopic studies of Cepheid variables at various phases similar 
to ours were already conducted by Andrievsky, Luck & Kov- 
tyukh (2005), Luck & Andrievsky (2004), Kovtyukh et al. 
(2005a), and Luck et al. (2008), for four different Cepheid 
groups of (i) 3 < P < 6, (ii) 6 < P < 10, (iii) 10 < P 
(where P is the pulsation period in day), and (iv) s-Cepheids 
with small amplitudes, respectively, in order to study the 
phase-dependent variations of stellar fundamental parame- 
ters. They showed that almost consistent abundances were 
obtained for various elements by using the spectroscopic pa- 
rameters determined based on Kovtyukh & Andrievsky's 
(2005) procedure. 

In contrast, we treat and discuss our targets (including 
classical Cepheids in a wide range of P as well as s-Cepheids) 
as a whole without classifying them into subgroups; but we 
particularly focus on determining the abundances of four 
elements (C, N, O, and Na) as precisely as possible by care- 
fully applying the spectrum-synthesis technique along with 
non-LTE corrections, which is the distinction compared to 
their studies. 

The remainder of this paper is organized as follows. We 
describe in Section 2 the observational data of the program 
stars (12 Cepheid variables). The spectroscopic determina- 
tion of atmospheric parameters corresponding to each spec- 
trum, which uses Fe I and Fe n lines, is explained in Sec- 
tion 3. In Section 4 are spelled out the procedures of our 
abundance analysis, which are made up of LTE spectrum- 
synthesis fitting, inverse evaluation of equivalent widths, and 
non-LTE abundance determinations from them. The results 
of atmospheric parameters as well as C, N, O, and Na abun- 
dances are discussed in Section 5 with respect to each of the 
specific check points enumerated above. The conclusion is 
summarized in Section 6. 



2 OBSERVATIONAL DATA 

Twelve representative Cepheid variables (including five s- 
Cepheids showing particularly small light-variation ampli- 
tude of | Ay | < 0.5 mag) were selected as our program stars, 
which are apparently bright (mostly V ~ 4-6 mag) and have 
a wide variety of pulsation periods (~2-16 days). The basic 
data of these targets are summarized in Table 1. 

The observations were carried out during the four nights 



on 2009 October 2-5 by using BOES (Bohyunsan Observa- 
tory Echelle Spectrograph) attached to the 1.8 m reflector at 
Bohyunsan Optical Astronomy Observatory. Using 2kx4k 
CCD (pixel size of 15 /mi x 15 |Um), this echelle spectro- 
graph enabled us to obtain spectra of wide wavelength cov- 
erage (from - 3800 A to - 9200 A). We used 200/mi fiber 
corresponding to the resolving power of R ~ 45000. The in- 
tegration time for each exposure was from a few minutes up 
to rsj 15-20 minutes depending on the brightness of a target. 
If ADU counts attained by one exposure were not sufficient, 
the second (or even third) exposure was tried for the same 
star and the successive frames were co-added, though only 
one exposure was enough in normal cases. In any case, the in- 
tegrated exposure time for each spectrum (whichever single 
or co-added) never exceeded 40 minutes, which guarantees 
that any serious blurring effect caused by line-profile varia- 
tions is unlikely. In such a way, each of the program stars 
were observed several times in a night with an interval of 
a few hours, though the actual frequency differed from star 
to star. Thus, as a result of 4-night observations, we could 
obtain a total of 122 spectra, which consist of 7-17 spectra 
per each star corresponding to different observational times. 

The reduction of the echelle spectra (bias subtraction, 
flat fielding, spectrum extraction, wavelength calibration, 
and continuum normalization) was carried out mainly with 
the software developed by Kang et al. (2006) and partly 
with IRAF^We could accomplish sufficiently high S/N ra- 
tio of several hundred at the relevant regions for most of 
the spectra (except for the last one of S Cep, for which the 
quality is considerably poor). Then, the apparent stellar ra- 
dial velocity was determined by comparing each spectrum 
in the orange region (A ~ 6100 A) with the template solar 
spectrum, which was further converted to the heliocentric 
radial velocity by applying the relevant correction evaluated 
with the help of IRAF task "rvcorrect." The fundamental 
information for each of the 122 spectra (observational time 
in Julian day, the corresponding pulsation phase <j>, and the 
heliocentric radial velocity V ra d) is presented in Table 2. 



3 ATMOSPHERIC PARAMETERS 

The atmospheric parameters necessary for constructing 
model atmospheres as well as for abundance determinations 
[T e ff (effective temperature), logg (surface gravity), £ (mi- 
croturbulent velocity dispersion), and [Fe/H] (metallicity, 
represented by the Fe abundance relative to the Sun)] were 
spectroscopically determined from the equivalent widths 
(W\) of Fe I and Fe n lines, based on the conventional 
requirements of (a) excitation equilibrium (Fe abundances 
show no systematic dependence on the excitation potential) , 

(b) ionization equilibrium (mean Fe abundance from Fe I 
lines and that from Fe II lines agree with each other), and 

(c) curve-of-growth matching (Fe abundances do not sys- 
tematically depend on W\). 

Practically, we used the program TGVIT (an updated 
version of the original program named TGV) developed for 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 
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this purpose (Takeda et al. 2005; cf. section 3.1 therein), 
which can establish these four parameters simultaneously 
by numerically finding the best solution such that minimiz- 
ing the dispersion function D 2 defined as the combination of 
(log ei) (mean Fe abundance derived from Fe I lines), (log €2) 
(mean Fe abundance derived from Fe 11 lines), o\ (standard 
deviation of the Fe abundances from Fe 1 lines) , and ai (stan- 
dard deviation of the Fe abundances from Fe 11 lines) as 

D 2 = (aj + a\) + «lo gei ) - <lo ge2 » 2 . (1) 

The details regarding the principle and algorithm of this 
method are described in Takeda, Ohkubo & Sadakane 
(2002). 

Regarding the observational W\ data of Fe lines, we 
could measure W\s of ~ 100-250 Fe I lines and ~ 15- 
25 Fe II lines on each of the 122 spectra (the number of 
measurable lines significantly depends on the spectral line 
width determined by macroscopic line broadening as well 
as on the S/N ratio), while consulting the line list (com- 
prising 302 Fe I lines and 28 Fe n lines) given in electronic 
table El of Takeda et al. (2005), where Gaussian fitting was 
adopted in most cases. In practice, avoiding the use of very 
strong lines showing appreciable damping wings, which are 
unsuitable for abundance determinations, we decided to use 
only those lines satisfying the condition of w < 200 mA, 
where w = W\ • (5000/ A) is the reduced equivalent width to 
A = 5000 A. 

We first carried out some test calculations by using the 
spectrum of 5 Cep (213306-091002A), in order to confirm the 
argument of Kovtyukh & Andrievsky (1999), who suggested 
that Fe n lines of various strengths can be safely used, while 
the mean abundance from Fe I lines should be determined 
based on only weak lines (as the weak- line limit). Namely, 
the following two different cases were examined concerning 
the input data set of W\ . 

— Case (A): 

The TGVIT program was applied without any specific con- 
straint on the line strengths (i.e., all available Fe I and Fe II 
lines satisfying w < 200 mA). Then, the solutions of [T e ff 
(K), logg (cm s _1 ), £ (km s _1 ), and [Fe/H] (dex)] minimiz- 
ing D 2 were obtained as [5694, 1.88, 3.56, and 0.00]. How- 
ever, the abundances from Fe I and Fe n lines corresponding 
to these parameters turned out appreciably W\ -dependent 
in an opposite manner (i.e., Fe I | Fe n abundances tend to 
progressively decrease | increase with an increase in W\) as 
shown in Fig. la, which is hardly acceptable. 

— Case (B): 

Next, we restricted Fe I lines to only those fairly weak ones 
satisfying the criterion of w ^ 30 mA, while Fe II lines are 
usually used as before, in analogy with treatment suggested 
by Kovtyukh & Andrievsky (1999). The resulting parame- 
ter solutions for this case were [5706, 2.19, 3.88, and +0.07]; 
the especially notable change compared to Case (A) is an 
increase of logg by ~ 0.3 dex. The loge vs. W\ relation for 
this case is displayed in Fig. lc, where we can satisfactorily 
confirm that the Fe abundances do not show any system- 
atic dependence upon W\. Naturally, as recognized in the 
loge vs. xiow plots, the abundance scatter around the mean 
is appreciably smaller for Case (B) (Fig. Id) in comparison 
with Case (A) (Fig. lb). 

Given these results, which reconfirmed the consequence 
of Kovtyukh & Andrievsky (1999), we decided to follow Case 



(B) and use only the weak Fe I lines with w ^ 30 mA but all 
Fe II lines with w < 200 mA, as the input W\ data to which 
TGVIT was applied. The number of actually used Fe I lines 
after this screening was ~ 40-130 (i.e., almost as half as 
the originally measured ones) . The finally resulting parame- 
ters (T e ff, log £, and [Fe/H]) for each of the 122 spectra are 
presented in Table 2. The typical statistical uncertainties in- 
volved in these solutions estimated in the manner described 
in Section 3.2 of Takeda et al. (2002) are (~ 50 K, - 0.1 dex, 
~ 0.3 km s _1 , and ~ 0.05 dex), respectively. The corre- 
sponding loge vs. W\ and loge vs. xw plots are depicted 
in Fig. 2 and Fig. 3, respectively, where we can see that the 
requirements (a), (b), and (c) are reasonably fulfilled. The 
detailed line- by- line data (W\, loge) for each stellar spec- 
trum are also given as the electronic data tables contained 
in the supplementary materials (the results for HD ?????? 
are presented in the file named "??????_feabw.dat"). 



4 ABUNDANCE DETERMINATIONS 

4.1 Target Lines of C, N, O, and Na 

Regarding the choice of lines to be used for deriving the 
abundances of C, N, O, and Na, several requirements were 
taken into consideration: 

— From the viewpoint of mutual consistency, lines of the 
same type had better be used for all these elements; so we 
would rely on permitted lines of neutral species. 

— The strength of the line must not be too weak, so 
that firm detectability in the relevant parameter ranges 
(T eff - 5000-7000 K, \ogg - 1-3) can be guaranteed. 

— On the other hand, we would like to avoid considerably 
strong lines, such as those very sensitive to a non-LTE effect 
or microturbulence. 

— Since our analysis is based on the spectrum synthesis 
method, it is preferable to select a wavelength region of 
moderate width comprising a few lines belonging the same 
multiplet, so that they may be analyzed at a time. 

We then decided to invoke the follow- 
ing lines which reasonably meet these condi- 
tions: C I 7111/7113/7115/7116/7119 lines for C, 
N I 8680/8683/8686 lines for N, O I 6155-8 lines for 
O, and Na I 6154/6161 lines for Na. 

4.2 Synthetic Spectrum Fitting 

Abundance determinations in the first step were carried 
out by using our synthetic spectrum fitting code MP- 
FIT, in which the spectrum synthesis part is originally 
based on Kurucz's (1993) WIDTH9 program. It accom- 
plishes the best-match between the theoretical and ob- 
servational spectra based on the algorithm described in 
Takeda (1995), by simultaneously varying the abundances 
of several key elements (logei, loge2, . . .), macrobroaden- 
ing parameter (vm), and the radial- velocity (wavelength) 
shift (AA). The macrobroadening parameter (vm) is the e- 
folding width of the Gaussian macrobroadening function, 
M(v) oc exp[— (v/vm) 2 ], which represents the combined ef- 
fects of instrumental broadening, macroturbulence, and ro- 
tational velocity (though it is essentially dominated by mac- 
troturbulence in the present case). 
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Practically, our spectrum fitting analysis was conducted 
for the following three wavelength regions, where the ele- 
ments whose abundances were treated as variables are enu- 
merated in each bracket: (i) 6143-6163 A [O, Na, Si, Ca, Fe] 
including O I 6155-8 and Na I 6154/6161 lines, (ii) 7110- 
7121 A [C, Fe, Ni] including C I 7111/7113/7115/7116/7119 
lines and (hi) 8677-8697 A [N, Si, S, Fe] including 
N I 8680/8683/8686 lines. 

Regarding the atomic data of spectral lines (wave- 
lengths, excitation potentials, oscillator strengths, and 
damping constants), we basically invoked the compilation 
of Kurucz & Bell (1995). However, pre- adjustments of sev- 
eral log gf values were necessary (i.e., use of empirically de- 
termined solar gf values) for several lines in 7110-7121 A 
region in order to accomplish a satisfactory match between 
the observed and theoretical spectrum. The finally adopted 
atomic data of important spectral lines are presented in ta- 
ble 2. The atmospheric model to be used for each spectrum 
was generated by interpolating Kurucz 's (1993) grid of AT- 
LAS9 model atmospheres in terms of T e ff, log and [Fe/H]. 

This fitting analysis turned out quite successful and we 
could obtain the abundances (especially for C, N, O, and 
Na) for all the 122 spectra. In case where some part of the 
spectrum was found to be damaged due to cosmic rays or 
telluric lines (e.g., occasionally encountered at A ~ 7118 A 
in the 7110-7121 A fitting), such a region was masked and 
excluded from the evaluation of ^{O — C) evaluation. How 
the theoretical spectrum corresponding to the converged so- 
lutions fits well with the observed spectrum is displayed in 
figure 4 (6143-6163 A fitting), figure 5 (7110-7121 A fitting), 
and figure 6 (8677-8697 A fitting). Note that we assumed 
LTE for all lines at this stage of synthetic spectrum- fitting. 
The solutions of the macrobroadening width (vm) derived 
from the 6143-6163 A fitting are given in Table 2, while fur- 
ther work is yet to be done to establish the final abundances 
as described in the next Section 4.3. 

4.3 Equivalent Widths Analysis 

Despite that the synthetic spectrum fitting directly yielded 
the abundance solutions of C, N, O, and Na (the main pur- 
pose of this study), this approach is not necessarily suit- 
able when one wants to examine how the results would be 
changed in different conditions (e.g., LTE vs. non-LTE, ap- 
plication of non-LTE calculations done at different assump- 
tions, abundance sensitivity to parameter changes, etc.), 
since it is tedious to repeat the fitting process to obtain the 
new solution. Therefore, with the help of the modified ver- 
sion of Kurucz's (1993) WIDTH9 prograrrj^] we computed 
the equivalent widths for C I lines (W7111, W7113, W7115, 
W7116, W7119), N 1 lines (Wseso, W 8 68 3 , Wsese), O 1 lines 
(HW- 8 )Q and Na 1 lines (W6154, W6i6i) "inversely" from 

2 This WIDTH9 program had been considerably modified by 
Y. Takeda in various respects; e.g., inclusion of non-LTE effects, 
treatment of total equivalent width for multi-component lines; 
etc. 

3 We use the total equivalent width of the O I 6155-8 feature 
consisting of 9 component lines (which appears as a merged triplet 
feature) in order to maintain consistency with Takeda & Takada- 
Hidai (1998) as well as Takeda et al. (2010). Since the actual O 1 
6155-8 feature is appreciably blended with lines of other elements 



the abundance solutions (resulting from spectrum synthe- 
sis) along with the adopted atmospheric model/parameters, 
which are much easier to handle. Based on such evaluated 
W\ values, the LTE (loge LTE ) as well as non-LTE abun- 
dances (loge NLTE ) were freshly computed, from which the 
non-LTE correction (A[= loge NLTE - loge LTE ]) was further 
derived. 

Regarding the calculations for evaluating the non-LTE 
departure coefficients, we followed the procedures described 
in Takeda & Takada-Hidai (2000) (for C), Takeda & Takada- 
Hidai (1995) (for N), Takeda & Takada-Hidai (1998) as well 
as Takeda et al. (2010) (for O), and Takeda & Takada-Hidai 
(1994) (for Na), which should be consulted for the details. 
In the actual non-LTE calculations for element X, two sets 
of departure coefficients were prepared corresponding to two 
different input theoretical abundances ([X/H]*, [X/H]^; ex- 
pressed as the values relative to the solar abundance) which 
were so chosen as to encompass the real abundances of the 
program stars; i.e., (-0.5, 0.0) for C, (0.0, +1.0) for N, (0.0, 
-0.5) for O, and (0.0, +0.5) for Na. Then, from a given 
W\, we obtained two kinds of non-LTE abundances (logei, 
log £2), which were further interpolated (or extrapolated) 
while requiring that the finally resulting non-LTE abun- 
dance be consistent with the input theoretical abundance 
(cf. Section 4.2 in Takeda & Takada-Hidai 1994). 

The values of the equivalent width (Wi), non-LTE abun- 
dance (logef LTE ), and non-LTE correction (A«) for each 
line i of C, N, O, and Na, which were derived for all 
the 122 spectra, are presented as the electronic data ta- 
bles (the results for HD ?????? are given in the file named 
«?????? _cnona.dat") contained in the supplementary mate- 
rials. 

The final non-LTE abundance and the non-LTE correc- 
tion for element X, loge(X) and Ax (X = C, N, O, and Na), 
were eventually obtained by averaging log eP™ 3 and A t over 
each relevant line i. We also calculated the differential abun- 
dance relative to the Sun as [X/H] = loge(X) — logeQ(X), 
which we will mainly use in the discussion. Regarding the 
reference solar abundances, we adopted the following val- 
ues, which were derived from the solar flux spectra based 
on practically the same lines and the same atomic data by 
taking into account the non-LTE effect: loge©(C) = 8.51 
(Takeda et al. 2013), loge (N) = 8.05 (newly derived for 
this study from the N 1 8683 line with the non-LTE correc- 
tion of —0.05 dex, for which the solar Ws683 was measured to 
be 6.1 mA), loge (O) = 8.81 (Takeda & Honda 2005), and 
log e (Na) = 6.32 (Takeda et al. 2003), where the abundance 
values are expressed in the usual normalization of log e(H) = 
12.00. The resulting [C/H], [N/H], [O/H], and [Na/H] (along 
with the corresponding Ac, An, Ao, and An & ) derived for 
each of the 122 spectra are summarized in Table 2. As seen 
in this table, the extents of the (negative) non-LTE correc- 
tions are not significant for C, O, and Na (< 0.1 dex), but 
rather considerable for N (~ 0.2-0.4 dex). 



(such as Si I, Ca I, and Fe 1 lines; cf. Fig. 4 and Table 3), this 
^6155—8) is n °t so much a realistic value measurable from actual 
spectra as rather an idealized quantity. 
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5 DISCUSSION 

5.1 Phase-Dependence of Physical Quantities 

The results of the physical parameters in the atmosphere 
(Vrad, T e ff, logg, £, vm) as well as the elemental abun- 
dances ([Fe/H], [C/H], [N/H], [0/H], and [Na/H]) are plot- 
ted against the phase {(f)) in Fig. 7-18 for each of the program 
stars (SU Cas, SZ Tau, RT Aur, £ Gem, FF Aql, 77 Aql, S Sge, 
X Cyg, T Vul, DT Cyg, V1334 Cyg, and 5 Cep). Although 
we do not discuss these figures for each of the individual stars 
in detail, we can recognize several typical characteristics al- 
ready known for Cepheid variables (see, e.g., Bersier, Burki 

6 Kurucz 1997), such as the enhancement of T e ff and turbu- 
lent velocity fields (£ and vm which are well- correlated) at 
the contraction or compressed phase of positive line-of-sight 
velocity (relative to the mean system velocity). 

In order to check the results of our analysis, our values 
of T e ff, logg, £, and [Fe/H] derived for S Cep (one of the 
best studied representative Cepheids) and those of previous 
studies (Luck & Lambert 1985; Kovtyukh, Komarov & De- 
penchuk 1994; Fry & Carney 1997; Kovtyukh & Andrievsky 
1999; Andrievsky, Luck & Kovtyukh 2005) are overplotted 
against <f> in Fig. 19, where we restricted the comparison to 
spectroscopically determined parameters derived in a simi- 
lar manner to ours. We can see several notable trends from 
this figure: 

— A remarkably good agreement with these literature val- 
ues is seen for T eff (Fig. 19a) and [Fe/H] (Fig. 19d). 

— Our spectroscopic gravities are almost consistent with 
Kovtyukh & Andrievsky 's (1999) "non-standard" results, re- 
flecting that both were derived in the same manner (mainly 
based on Fe 11 lines along with the limited use of only very 
weak Fe 1 lines ), while older results using stronger Fe 1 
lines (e.g., Fry & Carney 1997) tend to be appreciably lower 
by ~ 0.5 dex (Fig. 19b). This confirms the argument of 
Kovtyukh & Andrievsky (1999) that the spectroscopic logg 
based on the new "non-stadard" approach are preferable, 
since those old spectroscopic logg values of Cepheids (con- 
ventionally determined) were often found to be significantly 
lower than the dynamical values determined from the empir- 
ically estimated mass and radius (see, e.g., Luck & Lambert 
1985). 

— Regarding the microturbulence, our £ values tend to be 
systematically higher by ^ 1 km s _1 in comparison with 
other literature results (Fig. 19c). The reason for this dis- 
crepancy is not clear, since it is observed not only in com- 
parison with older work (where Fe 1 lines were used to deter- 
mine £) but also with the non-standard results of Kovtyukh 
& Andrievsky (1999) (where £ was determined by Fe 11 lines 
as in this study). 



5.2 Correlation with Pulsation Period 

Given that our program stars were so chosen as to cover 
a wide range of pulsation periods (P ~ 2-16 days), it is 
worth examining how our spectroscopically determined at- 
mospheric parameters depend on P, which is closely cor- 
related with the absolute magnitude (also with the stellar 
mass) through the period-luminosity relation. Fig. 20 shows 
all the results of T e ff , log g, £, vm, and [Fe/H] derived for each 
of the 12 stars, plotted against log P. 



We can see a clear trend of T e ff decreasing with an in- 
crease in logP (Fig. 20a), which can be roughly represented 
by the linear relation 

log T efT = 3.84 - 0.10 log P (2) 

(T e ff is in K and P is in days) with an uncertainty of a few 
hundred K. This reflects the fact the Cepheid instability 
strip is tilted in the theoretical logL vs. logT e ff diagram, 
in the sense that higher- L Cepheids with longer P tend to 
have lower T e ff(^] 

We then discuss the P-dependence of logg (Fig. 20b). 
The surface gravity is expressed in terms of M, L, and T e ff 
as 

logte/00) = log(M/M )-log(L/L )+41og(Teff/Teff, ).(3) 

If we use Benedict et al.'s (2007) period-luminosity relation 
for galactic Cepheids (cf. their fig. 5) 

M v = -1.62 - 2.43 log P (4) 
we have 

log(L/L©) = 2.54 - 0.4P.C. + 0.97 log P, (5) 

where B.C. is the bolometric correction. We then invoke Ca- 
puto et al.'s (2005) mass-luminosity relation for solar com- 
position (cf. their eq. (2)) 

log(L/L ) = 0.72 + 3.35 log(M/M ). (6) 

Combining Eq. (5) and (6), we get a mass-period relation 
as 

log(M/M ) = 0.54 - 0.12P.C. + 0.29 log P. (7) 

Consequently, inserting Eq. (2), (5), and (7) into Eq. (3), we 
finally obtain 

log g = 2.75 + 0.28P.C. - 1.08 log P ~ 2.72 - 1.08 log P (8) 

as the P-dependence of dynamical logg, where we assumed 
B.C. ~ —0.1 (actual values range from ~ —0.2 to ~ 0.0 
depending on T e ff, but insignificant; e.g., Flower 1996). This 
relation is also depicted in Fig. 20b. As seen from this figure, 
the spectroscopic logg values derived in this study tend to 
be systematically larger by this relation of dynamical logg, 
which indicates that some of the relations we assumed above 
may not be appropriate or some systematic error may be in- 
volved in our logg. Nevertheless, considering the large dis- 
persion of observed logg along with the incomplete phase 
coverage of our data, we may regard that both are in tolera- 
ble consistency, despite the discrepancy especially in short- 
period Cepheids of P ^ 2-3 days. 

Fig. 20c and 20d suggest that £ as well as vm attain 
an apparent maximum at logP ^0.8 and turn to decrease 
with increasing P, which may contradict the intuitive expec- 
tation that the atmospheric turbulence would progressively 
grow with a decrease of log g because of increased instability 

4 Although Luck et al. (2008) concluded that the mean T e ff for s- 
Cepheids tends to be higher than that of classical Cepheids, such 
a trend can not be confirmed in Fig. 20a for the five s-Cepheids 
(SU Cas, SZ Tau, FF Aql, DT Cyg, and V1334 Cyg) included in 
our targets, which might be due to the incomplete phase coverage 
as well as the insufficient number of our program stars compared 
with their study (9 s-Cepheids and 30 classical Cepheids). 
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in lower- density condition. We consider, however, that this 
should not be taken too seriously, since the phase coverage 
of X Cyg (having the longest P of 16 days) is insufficient and 
our 9 spectra for this star correspond to the expanding and 
decelerating phase (<j) 0.2-0.4) of non-enhanced turbulence 
(Fig. 14). 

The dispersion of [Fe/H] for each star in Fig. 20e sug- 
gests that the precision of its determination for a spectrum 
at a given phase (as a by-product of atmospheric parameter 
determination based on Fe I and Fe n lines) is on the order of 
~ 0.1 dex, though with rather large star-to-star differences 
(e.g., the reason for the appreciably larger [Fe/H] disper- 
sion in S Sge is because our data correspond to <j) ~ 0.7- 
1.1, where the spectral variation is rapid and the analysis 
is comparatively more difficult). Overall, these [Fe/H] val- 
ues scatter around ~ without any systematic dependence 
of [Fe/H] upon P, which means that all the program stars 
have almost the solar metallicity irrespective of the absolute 
magnitude and the mass. 

Somewhat unexpectedly, such P-independent chemical 
homogeneity (i.e., small star-to-star abundance dispersion) 
was found also for the other elements. By averaging the 
[X/H]i results for each individual spectrum z, we computed 
the mean abundance ([X/H]) and the standard deviation 
ax for each of the 12 stars (X = C, N, O, Na, and Fe), 
which are given in Table 2 (expressed in italic at the first 
line of each star's section). These ([C/H]), ([N/H]), ([O/H]), 
([Na/H]), and ([Fe/H]) for each star are plotted against log P 
in Fig. 21a-e, where the extents of relevant ±a are shown 
as error bars. These figures demonstrate that none of these 
([X/H])s show any significant P-dependence, which are prac- 
tically homogeneous within a dispersion of < 0.1-0.2 dex. 
We may thus conclude ([C/H]) - -0.3, ([N/H]) - +0.3- 
0.4, ([O/H]) - 0, ([Na/H]) - +0.2, and ([Fe/H]) - as the 
abundance characteristics for these elements, which hold for 
all the program stars regardless of log P. We note that these 
values are almost consistent with the recent extensive stud- 
ies for various Cepheid variables (Andrievsky et al. 2005; 
Luck & Andrievsky 2004; Kovtyukh et al. 2005a; Luck et al. 
2008), where we can see that their C, N, O, and Na abun- 
dances are almost the P-independent with the means (and 
the standard deviations) of ([C/H]) = -0.15 (a = 0.08; 
38 stars excluding SV Mon showing an exceptionally large 
deficiency of -0.90), ([N/H]) = +0.38 (a = 0.17; 6 stars), 
([O/H]) = -0.08 (a = 0.09; 38 stars), and ([Na/H]) = +0.21 
(a = 0.07; 39 stars). 

5.3 Abundance Trends and Their Implications 

Now, since the abundances of key elements have been estab- 
lished, we can discuss the main problem which motivated 
this investigation (cf. Section 1): the nature of envelope mix- 
ing inferred from the characteristics of C, N, O, and Na 
abundances. 

According to the results derived in Section 5.2, C is 
mildly deficient by ~ 0.3 dex, N is appreciably enhanced by 
~ 0.3-0.4 dex, Na is moderately overabundant by ~ 0.2 dex, 
while O (as well as Fe) is practically solar without any signif- 
icant peculiarity (within an uncertainty of ~ 0.1 dex). This 
tendency can be visually confirmed in Fig. 22a and Fig. 22b, 
where ([N/H]) vs. ([C/H]) and ([Na/H]) vs. ([O/H]) corre- 
lations are plotted, respectively. 



Regarding the oxygen problem to be clarified in the first 
place, we can confidently state that any significant mixing 
of ON-cycle product has not taken place in the envelope of 
Cepheid variables. This means that the result of apparently 
subsolar [O/H] derived by Luck & Lambert (1985) is most 
likely due to the use of inappropriate atmospheric param- 
eters (especially too low logg), as suggested by Kovtyukh 
& Andrievsky (1999). Then, the abundance peculiarities of 
C and N must have been caused mainly by the dredge-up 
of CN-cycle product. The weak ant i- correlation between the 
C and N abundances implied from Fig. 22a is qualitatively 
consistent with this scenario. 

In order to check this quantitatively, the CNO abun- 
dances (loge) of our program stars are plotted in the 
(loge(N)) vs. (loge(C+0)) diagram and the (loge(N)) vs. 
(loge(C)) diagram in Fig. 22c and Fig. 22d, respectively, 
where the expected relations under the condition that 
e(C+N+0) is conserved are also drawn. Further, the re- 
sulting sums of loge(C+N+0) evaluated for the program 
stars are plotted against logP in Fig. 22e. We can confirm 
that the observed abundances are almost on the expected 
curve in Fig. 22c (compare this figure with Fig. 1 of Luck & 
Lambert 1985), and that loge(C+N+0) is almost conserved 
at the solar value (Fig. 22f). 

It should be remarked, however, that we can not rule 
out the possibility of only a slight underabundance in O 
caused by the dredge-up of some (not significant) ON-cycle 
product. Actually, Fig. 22d implies that loge(N) vs. loge(C) 
correlation can not necessarily be well described only by the 
conservation conditions of loge(C+N) = loge©(C+N) and 
loge(O) = log e (O) for stars of loge(N)> 8.5. Rather, it 
is more reasonable to consider that a mixing of some ON- 
cycled material could have caused a slight O-deficiency by 
several hundredths dex (within < 0.1 dex) in order to ex- 
plain the loge(N) values of these stars (Fig. 22d). 

Finally, our results of moderate enhancement of Na by 
^0.2 dex indicate that the NeNa-cycle product is dredged- 
up in the envelope of Cepheid variables. Takeda & Takada- 
Hidai (1994) derived [Na/H] = +0.15 (rj Aql) and [Na/H] 
= —0.07 (C Gem) based on their non-LTE analysis of the 
Na I 8194.12 line. Comparing these with the present results 
of +0.16 (77 Aql) and +0.23 (£ Gem), we can see an appre- 
ciable discrepancy for the latter. However, their old analysis 
should be regarded as less reliable, where the atmospheric 
parameters corresponding to the observed phase were not 
properly determined but roughly assumed while consulting 
the literature values. 

5.4 Mass-Independent Mixing in Cepheids 

Before starting this investigation, we anticipated that a sig- 
nificantly large abundance scatter may be found at least for 
C and N (or possibly also for O), as reported in the pio- 
neering work by Luck & Lambert (1985; cf. Fig. 1 or Fig. 9 
therein), and that such a diversity (if exists) might be corre- 
lated with P (or equivalently, the stellar mass M). However, 
the considerably small star-to-star dispersion in the abun- 
dances of all elements (C, N, O, Na, Fe) was rather an unex- 
pected result. Though, admittedly, a weak anti-correlation 
between N and C abundances (Fig. 22a, Fig. 22d) suggests 
that the extent of dredged-up material mixed in envelope 
may differ slightly from star to star, it does not appear to 
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have anything to do with M, since [N/C] does not show any 
systematic dependence upon P (Fig. 22e). So, another im- 
portant outcome of this study is the confirmation that the 
extent of mixing in the envelope of Cepheid variables does 
not depend upon the stellar mass. Interestingly, essentially 
the same conclusion was reached by Kovtyukh, Wallerstein 
& Andrievsky (2005b), who carried out an extensive spec- 
troscopic abundance study for 16 distant Cepheids with a 
wide range of P (~ 3-27 days) and found that neither [C/N] 
(~ +0.6) nor [Na/Fe] (~ +0.2) show any systematic depen- 
dence upon logP, which is quite consistent with our results. 

We should note, however, that this consequence is guar- 
anteed only for the pulsation variables within the Cepheid 
instability strip, which should not be simply extended to 
non-variable supergiants or giants in different parameter 
ranges (e.g., in terms of T e ff or M). For example, we can 
not say much about the argument of Takeda & Takada- 
Hidai (1994) (Na tends to be more enriched with an in- 
crease of mass for supergiants of ~ 10-25 M©; cf. Fig. 7 
therein), since the mass range of the program stars in this 
study is between ~ 4M© and ~ 8M© according to Eq. (7), 
and their sample was mostly non-variable supergiants. Sim- 
ilarly, the systematic log ^-dependence of [Na/H] concluded 
by Andrievsky et al. (2002) was mainly due to supergiants 
(not Cepheids) of 0.5 < logg < 2 (cf. their fig. 1), which 
is again less relevant to this study. Besides, Smiljanic et al. 
(2006) concluded based on their CNO abundance study for 
19 no n- variable supergiants covering a rather large range of 
stellar parameters (4100 K < T eff < 7500 K, 0.8 < \ogg < 
2.5, 2 < M/Mq < 13) that the [N/C] ratios tend to increase 
with M and have considerably larger values (up to < +2) 
than the theoretical expectation. Further, Takeda, Sato & 
Murata (2008) reported in their analysis of 322 late-G and 
early-K giants (4500 K < T eff < 5500 K, 1.5 < \ogg < 3.5, 1 
< M/Mq < 5) that [C/Fe], [O/Fe], and [Na/Fe] show subso- 
lar, subsolar, and supersolar tendency, respectively, with the 
degree of peculiarity increasing with M (cf. Fig. 12 therein), 
which may imply that a mass-dependent dredge-up of CN- 
cycle, ON-cycle, and NeNa-cycle products may take place in 
the envelope of red giants. 

Consequently, we would still consider it quite possi- 
ble that the mixing-induced abundance peculiarities of non- 
variable intermediate-mass evolved stars in general are sig- 
nificantly diversified and may depend on M. If so, we would 
have to find a reasonable answer to the question "Why the 
absence of large star-to-star dispersion or of M-dependence 
in the evolution-induced mixing is limited only to Cepheid 
variables?" , to which contributions by theoreticians are de- 
sirably awaited. 



photospheric abundances of C, N, O, and Na, which are the 
key elements for investigating the dredge-up of H-burning 
products from the interior, based on 122 high-dispersion 
spectra (~ 10 spectra of different phases per target) of wide 
wavelength coverage collected at Bohyunsan Astronomical 
Observatory. 

The atmospheric parameters (T e ff, logg, £, and [Fe/H]) 
corresponding to each phase were determined spectroscopi- 
cally from the equivalent widths of Fe I and Fe n lines by the 
requirements of excitation equilibrium, ionization equilib- 
rium, and curve-of- growth matching. Practically, we applied 
Takeda et al.'s (2005) TGVIT program, while using only 
fairly weak Fe ilines (< 30 mA) as suggested by Kovtyukh & 
Andrievsky (1999) to avoid non-LTE-sensitive stronger Fe I 
lines. The resulting parameters as well as the radial veloci- 
ties were confirmed to show the well-known phase-dependent 
characteristics of Cepheids. 

The abundances of C, N, O, and Na were derived by 
applying the spectrum-synthesis fitting technique to three 
wavelength regions (6143-6163 A, 7110-7121 A, and 8677- 
8697 A) comprising C I 7111/7113/7115/7116/7119, O I 
6155-8, N I 8680/8683/8686, and Na I 6154/6161 lines. 
Then, from the equivalent widths inversely computed from 
the fitting-based abundances, the final abundances were ob- 
tained by taking into account the non-LTE effect. 

The resulting abundances of these elements for 12 pro- 
gram stars turned out to show remarkably small star-to-star 
dispersions (< 0.1-0.2 dex) without any significant depen- 
dence upon the pulsation period: near-solar Fe ([Fe/H] ~ 
0.0), moderately underabundant C ([C/H] ~ —0.3), appre- 
ciably overabundant N ([N/H] ~ +0.4-0.5), and mildly su- 
persolar Na ([Na/H] ~ +0.2). Based on these results, we 
conclude as follows. 

— (1) These CNO abundance trends can be interpreted 
mainly as due to the canonical dredge-up of CN-cycled mate- 
rial, while the significant non-canonical deep mixing of ON- 
cycled gas is ruled out (though only a slight mixing may still 
be possible). 

— (2) The mild but definite overabundance of Na suggests 
that the NeNa-cycle product is also dredged up. 

— (3) The extent of mixing-induced peculiarities in the en- 
velope of Cepheid variables is almost independent on My as 
well as on M. However, given the observational suggestion 
that a significant diversity or a tendency of M-dependence 
exists in the C, N, O, and Na abundances of no n- variable 
supergiants/giants of other types, the problem of "why such 
a homogeneity in the surface abundances is limited only to 
Cepheids" is yet to be investigated further. 



6 CONCLUSION 

The mixing in the envelope of intermediate-mass F-G su- 
pergiants causing surface abundance peculiarities by the 
dredge-up of H-burning products is not yet well understood. 
Especially, regarding oxygen, whether or not its anomaly ex- 
ists due to non-canonical dredge-up of ON-cycled material 
is still controversial. 

In order to clarify this situation, we conducted an exten- 
sive spectroscopic study for selected 12 Cepheid variables of 
various pulsation periods (~ 2-16 days), and determined the 
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Table 1. Basic data of the program stars. 



Name 


HD# 


«2000 


#2000 


Cep.Type 


Sp.Type 


V 


Epoch 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SU Cas 


017463 


02:51:58.8 


+68:53:19 


s 


F5Ib/II-F7Ib/II 


5.70-6.18 


50100.156 


1.949319 


SZ Tau 


029260 


04:37:14.8 


+18:32:35 


s 


F5Ib-F9.5Ib 


6.33-6.75 


50101.605 


3.14873 


RT Aur 


045412 


06:28:34.1 


+30:29:35 


cl 


F41b-Gllb 


5.00-5.82 


50101.159 


3.728115 


C Gem 


052973 


07:04:06.5 


+20:34:13 


cl 


F7Ib-G3Ib 


3.62-4.18 


50108.93 


10.15073 


FF Aql 


176155 


18:58:14.7 


+17:21:39 


s 


F5Ia-F8Ia 


5.18-5.68 


50102.387 


4.470916 


r\ Aql 


187929 


19:52:28.4 


+01:00:20 


cl 


F6Ib-G4Ib 


3.48-4.39 


50100.861 


7.176641 


S Sge 


188727 


19:56:01.3 


+16:38:05 


cl 


F6Ib-G5Ib 


5.24-6.04 


50105.348 


8.382086 


XCyg 


197572 


20:43:24.2 


+35:35:16 


cl 


F7Ib-G8Ib 


5.85-6.91 


50106.014 


16.386332 


T Vul 


198726 


20:51:28.2 


+28:15:02 


cl 


F5Ib-G0Ib 


5.41-6.09 


50101.410 


4.435462 


DT Cyg 


201078 


21:06:30.2 


+31:11:05 


s 


F5.5-F7Ib/II 


5.57-5.96 


50102.487 


2.499215 


V1334 Cyg 


203156 


21:19:22.2 


+38:14:15 


s 


F2Ib 


5.77-5.96 


50102.549 


3.332816 


S Cep 


213306 


22:29:10.3 


+58:24:55 


cl 


F51b-Gllb 


3.48-4.37 


50102.860 


5.366341 



Following the star name and HD number (serial number in the Henry-Draper catalog) in Columns 
(1) and (2), the data of star coordinates (right ascension in HH:MM:SS.S and declination in 
deg:arcmin:arcsec), Cepheid type ("cl" for classical Cepheids and "s" for s-Cepheids), spectral type, 
and apparent visual magnitude (in mag) are presented in Columns (3)-(7), which were taken from 
the web site of General Catalogue of Variable Stars (http://www.sai.msu.su/gcvs/gcvs/index.htm). 
The last two Columns (8) and (9) give the epoch of brightness maximum expressed in Julian day 
(JD —2400000) and the pulsation period (in day), for which we consulted Table 1 of Kiss (1998). 
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Table 2. Results of stellar parameters and abundances derived from each spectrum of different phase. 



Code JD V rad T eff logs [Fe] £ v M [C] A c [N] A N [O] A Q [Na] A Na 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) 



SU Cas (HD 017463) 
















+ 0.06(±0.09) 






-0.30(±0. 


05) 


+ 0.46(±0. 


05) 


+ 


.05(±0. 


06) 


+ 0.20(±0. 


03) 


017463-091002A 


6.99 


0. 


.503 


— 1 


,4 


6500 


3. 


.43 


+o. 


.27 


3.6 


7. 


.6 


-0.26 


-0, 


.04 


+0.55 


-0. 


.26 


+o. 


.11 


-0. 


.03 


+0. 


28 


-0 


.07 


017463-091002B 


7.22 


0. 


.623 


+ !• 


3 


6336 


3. 


.04 


+o. 


.14 


3.7 


9. 


.1 


-0.27 


-0 


.05 


+0.48 


-0. 


,27 


+o. 


.09 


-0 


.03 


+ 0. 


21 


-0 


.07 


017463-091002C 


7.30 


0. 


.663 


+ 1. 


,3 


6469 


3. 


46 


+o. 


.22 


4.2 


9. 


.6 


-0.23 


-0. 


.04 


+0.55 


-0. 


.24 


+o. 


.15 


-0. 


,03 


+o. 


25 


-0 


.07 


017463-091003A 


7.98 


0. 


.012 


— 14, 


,8 


6833 


2. 


.76 


-0. 


.08 


3.9 


10. 


.9 


-0.38 


-0. 


.07 


+0.38 


-0 


,41 


— 0. 


.03 


-0. 


.05 


+o. 


16 


-0 


.07 


017463-091003B 


8.08 


0. 


.062 


— 15, 


.8 


6935 


2. 


.88 


-0 


.01 


3.8 


10. 


.9 


-0.37 


-0. 


.07 


+0.36 


-0 


.41 


-0. 


.07 


-0, 


.05 


+o. 


20 


-0 


.07 


017463-091003C 


8.16 


0. 


.103 


— 14, 


.9 


6853 


2. 


.89 


+o. 


.00 


3.7 


10. 


,4 


-0.36 


-0. 


.07 


+0.40 


-0. 


,41 


-0, 


.04 


-0 


.05 


+o. 


21 


-0 


.07 


017463-091003D 


8.21 


0. 


.132 


— 15. 


.0 


6916 


3. 


.19 


+o. 


.02 


4.1 


9. 


.7 


-0.31 


-0 


.06 


+0.43 


-0. 


,36 


+o. 


.01 


-0 


,04 


+o. 


23 


-0 


.07 


017463-091003E 


8.26 


0. 


.156 


— 14. 


,0 


6632 


2. 


,75 


+o. 


.00 


3.5 


9. 


.6 


-0.33 


-0. 


.07 


+0.45 


-0. 


.39 


+o. 


05 


-0. 


,05 


+o. 


16 


-0 


.07 


017463-091003F 


8.30 


0. 


.178 


— 13, 


,6 


6599 


2. 


.75 


+o. 


.03 


3.5 


9. 


.3 


-0.33 


-0. 


.07 


+0.45 


-0. 


.39 


+o. 


.03 


-0 


.05 


+o. 


17 


-0 


.07 


017463-091003G 


8.34 


0. 


.197 


— 12 


,7 


6660 


2. 


.86 


+o. 


.04 


3.5 


9. 





-0.35 


-0 


.07 


+0.44 


-0 


.38 


+o. 


.02 


-0 


,05 


+o. 


21 


-0 


.07 


017463-091004A 


9.05 


0. 


.563 


+o. 


,3 


6279 


2. 


.83 


+0. 


.08 


3.5 


7. 


5 


-0.27 


-0. 


,06 


+0.46 


-0. 


.30 


+o. 


05 


-0. 


.04 


+o. 


19 


-0 


.07 


017463-091004B 


9.14 


0. 


.605 


+ 1. 


,7 


6279 


2. 


,90 


+o. 


.09 


3.6 


8. 


,0 


-0.24 


-0 


.06 


+0.47 


-0. 


,28 


+o. 


.06 


-0. 


.04 


+o. 


19 


-0 


.07 


017463-091004C 


9.21 


0. 


.643 


+ 2. 


,6 


6278 


2. 


.86 


+o. 


.08 


3.7 


8. 


5 


-0.28 


-0. 


.06 


+0.47 


-0. 


.29 


+o. 


06 


-0. 


.04 


+o. 


.18 


-0 


.07 


017463-091004D 


9.28 


0. 


.677 


+ 2. 


.6 


6410 


3. 


.29 


+o. 


.17 


4.1 


9. 


.0 


-0.20 


-0. 


.04 


+0.54 


-0 


.25 


+o. 


.12 


-0 


.03 


+o. 


23 


-0 


.07 


017463-091005A 


9.97 


0. 


.034 


-15, 


.0 


6964 


3. 


.14 


+o. 


.01 


4.1 


9. 


.5 


-0.32 


-0. 


.06 


+0.42 


-0 


.38 


+o. 


.01 


-0 


.04 


+o. 


19 


-0 


.07 


017463-091005B 


10.22 


0. 


.160 


-14. 


,6 


6596 


2. 


,78 


+o. 


.00 


3.8 


8. 


.8 


-0.30 


-0. 


.07 


+0.47 


-0. 


,38 


+o. 


09 


-0. 


,05 


+o. 


16 


-0 


.07 


017463-091005C 


10.28 


0. 


.190 


-13. 


.3 


6629 


2. 


.90 


+o. 


.04 


3.6 


8. 


.5 


-0.30 


-0. 


.07 


+0.47 


-0 


,37 


+o. 


.08 


-0. 


,05 


+o. 


19 


-0 


.07 


SZ Tau (HD 029260) 
















+o 


.05(±0.04) 






-0.26(±0. 


04) 


+ 0.48(±0. 


05) 


+o 


.03(±0. 


04) 


+o 


.22(±0. 


03) 


029260-091002A 


7.15 


0. 


.704 


+ 10. 


.0 


6010 


2. 


,42 


+o. 


.01 


4.3 


12. 


,2 


-0.24 


-0 


.07 


+0.54 


-0, 


.31 


+o. 


,07 


-0. 


.05 


+o. 


18 


-0 


.08 


029260-091002B 


7.24 


0. 


.732 


+ 8. 


.9 


6143 


2. 


.71 


+o. 


.07 


4.9 


13. 


,0 


-0.22 


-0. 


.06 


+0.57 


-0 


.28 


+o. 


.11 


-0 


.04 


+o. 


22 


-0 


.07 


029260-091002C 


7.32 


0. 


.758 


+7. 


.2 


6141 


2. 


,42 


-0. 


.01 


4.5 


13. 


,7 


-0.29 


-0. 


.07 


+0.46 


-0. 


,32 


-0. 


,01 


-0. 


,05 


+o. 


18 


-0 


.07 


029260-091003A 


8.17 


0. 


.027 


-11. 


.0 


6405 


2. 


.41 


+o. 


.06 


4.0 


12. 


.1 


-0.31 


-0. 


.08 


+0.41 


-0. 


.38 


-0 


,03 


-0. 


,06 


+o. 


25 


-0 


.07 


029260-091003B 


8.23 


0. 


.045 


-10, 


.6 


6272 


2. 


,24 


+o. 


.00 


3.9 


11. 


.6 


-0.29 


-0. 


.09 


+0.46 


-0 


.39 


+o. 


.01 


-0 


.07 


+o. 


21 


-0 


.08 


029260-091003C 


8.27 


0. 


.059 


-9 


.7 


6276 


2. 


.35 


+o. 


.01 


4.0 


11. 


,4 


-0.25 


-0. 


.09 


+0.48 


-0 


.38 


+o. 


.04 


-0 


.06 


+o. 


21 


-0 


.08 


029260-091003D 


8.32 


0. 


.073 


-9 


.8 


6427 


2. 


.81 


+o. 


.12 


4.5 


11. 


,2 


-0.17 


-0. 


.06 


+0.53 


-0 


.33 


+o. 


.07 


-0 


.05 


+o. 


28 


-0 


.07 


029260-091004A 


9.15 


0. 


.337 


-0. 


.8 


6044 


2. 


,50 


+o. 


.08 


3.9 


7. 


.7 


-0.25 


-0. 


.07 


+0.45 


-0. 


,29 


-0. 


,04 


-0. 


,04 


+o. 


25 


-0 


.08 


029260-091004B 


9.22 


0. 


.362 


+o. 


.6 


5968 


2. 


.46 


+o. 


.06 


3.8 


7. 


.6 


-0.20 


-0. 


.07 


+0.47 


-0 


.28 


+o. 


02 


-0 


.04 


+o. 


23 


-0 


.08 


029260-091004C 


9.29 


0. 


.384 


+1. 


.9 


5919 


2. 


.30 


+o. 


.03 


3.8 


7. 


.7 


-0.26 


-0, 


.07 


+0.47 


-0. 


.29 


-0. 


.01 


-0 


.05 


+o. 


21 


-0 


.08 


029260-091005A 


10.15 


0. 


.655 


+11. 


.7 


5961 


2. 


.35 


+o. 


.06 


3.9 


11. 


.1 


-0.26 


-0 


.07 


+0.46 


-0 


.30 


+o. 


.01 


-0 


.05 


+o. 


21 


-0 


.08 


029260-091005B 


10.30 


0. 


.704 


+ 10. 


.4 


6063 


2. 


.36 


+o. 


.06 


4.3 


12. 


,2 


-0.32 


-0. 


.08 


+0.42 


-0 


.31 


+o. 


.05 


-0 


.05 


+o. 


21 


-0 


.08 



RT Aur (HD 045412) 










+ 0.09(±0.07) 




-0.25(±0.07) 


+ 0.47(±0.05) 


-0.01(±0.05) 


+ 0.25(±0.05) 


045412-091002A 


7.18 


0.776 


+37.2 


5695 


2.08 


+0.01 


4.6 


10.5 


-0.20 


-0.07 


+0.44 


-0.25 


+0.06 


-0.04 


+0.18 


-0.07 


045412-091002B 


7.27 


0.800 


+37.1 


5746 


2.18 


+0.08 


4.6 


10.8 


-0.23 


-0.07 


+0.41 


-0.24 


-0.01 


-0.04 


+0.20 


-0.07 


045412-091003A 


8.18 


0.043 


+3.3 


7001 


2.81 


+0.06 


4.5 


10.1 


-0.34 


-0.07 


+0.45 


-0.45 


-0.08 


-0.05 


+0.23 


-0.06 


045412-091003B 


8.24 


0.060 


+ 2.8 


7000 


2.85 


+0.02 


5.0 


10.0 


-0.35 


-0.07 


+0.45 


-0.43 


-0.06 


-0.05 


+0.22 


-0.06 


045412-091003C 


8.29 


0.072 


+ 1.7 


6890 


2.61 


+0.00 


4.4 


10.3 


-0.35 


-0.08 


+0.46 


-0.46 


-0.05 


-0.05 


+0.22 


-0.06 


045412-091003D 


8.33 


0.083 


+ 1.6 


7005 


3.03 


+0.10 


4.8 


10.1 


-0.27 


-0.06 


+0.51 


-0.42 


-0.03 


-0.04 


+0.24 


-0.06 


045412-091004A 


9.16 


0.306 


+ 12.3 


6269 


2.75 


+0.17 


3.5 


6.7 


-0.19 


-0.07 


+0.54 


-0.32 


+0.05 


-0.04 


+0.31 


-0.08 


045412-091004B 


9.25 


0.329 


+ 13.6 


6223 


2.69 


+0.18 


3.5 


6.2 


-0.20 


-0.07 


+0.51 


-0.32 


+0.01 


-0.04 


+0.30 


-0.08 


045412-091004C 


9.32 


0.349 


+ 15.0 


6250 


2.92 


+0.22 


3.6 


5.8 


-0.16 


-0.06 


+0.55 


-0.29 


+0.04 


-0.04 


+0.32 


-0.07 


045412-091005A 


10.24 


0.595 


+ 28.8 


5843 


2.43 


+0.10 


4.0 


7.4 


-0.21 


-0.06 


+0.42 


-0.24 


-0.01 


-0.04 


+0.24 


-0.07 


C Gem (HD 052973) 












+0.04 (±0.03) 




-0.29(±0.10) 


+ 0.35(±0.05) 


-0.14(±0.11) 


+ 0.23 


(±0.03) 


052973-091002A 


7.20 


0.405 


+ 17.5 


5321 


1.68 


+0.05 


4.4 


8.2 


-0.17 


-0.07 


+0.38 


-0.19 


-0.14 


-0.03 


+0.19 


-0.07 


052973-091002B 


7.29 


0.413 


+ 18.3 


5342 


1.90 


+0.06 


4.4 


8.4 


-0.12 


-0.06 


+0.40 


-0.18 


+0.09 


-0.03 


+0.19 


-0.07 


052973-091003A 


8.25 


0.508 


+ 22.6 


5452 


1.76 


+0.02 


4.7 


11.2 


-0.28 


-0.07 


+0.29 


-0.20 


-0.13 


-0.04 


+0.24 


-0.07 


052973-091003B 


8.29 


0.512 


+ 22.6 


5462 


1.72 


+0.04 


4.6 


11.3 


-0.29 


-0.08 


+0.29 


-0.21 


-0.14 


-0.04 


+0.25 


-0.07 


052973-091003C 


8.33 


0.516 


+ 22.5 


5501 


1.78 


+0.06 


4.5 


11.4 


-0.27 


-0.08 


+0.28 


-0.22 


-0.22 


-0.04 


+0.26 


-0.07 


052973-091004A 


9.25 


0.607 


+ 17.5 


5691 


2.08 


+0.09 


5.1 


11.9 


-0.33 


-0.07 


+0.38 


-0.23 


-0.20 


-0.04 


+0.27 


-0.07 


052973-091004B 


9.33 


0.615 


+ 16.9 


5581 


1.62 


-0.01 


4.6 


11.8 


-0.41 


-0.09 


+0.35 


-0.26 


-0.32 


-0.05 


+0.22 


-0.07 


052973-091005A 


10.25 


0.706 


+5.8 


5688 


1.79 


+0.00 


4.8 


9.8 


-0.41 


-0.09 


+0.43 


-0.28 


-0.07 


-0.05 


+0.21 


-0.07 


FF Aql (HD 176155) 










+ 0.04 (±0.06) 




-0.34(±0.03) 


+ 0.54(±0.02) 


-0.06(±0.02) 


+ 0.27(±0.02) 


176155-091002A 


6.93 


0.356 


-18.8 


6246 


2.57 


+0.16 


5.1 


8.0 


-0.27 


-0.07 


+0.58 


-0.33 


-0.02 


-0.05 


+0.31 


-0.07 


176155-091003A 


7.92 


0.578 


-9.7 


6081 


2.28 


+0.02 


5.5 


9.1 


-0.35 


-0.08 


+0.52 


-0.32 


-0.07 


-0.05 


+0.25 


-0.07 


176155-091003B 


7.99 


0.593 


-9.7 


6099 


2.29 


+0.03 


5.5 


9.2 


-0.34 


-0.08 


+0.51 


-0.32 


-0.09 


-0.05 


+0.26 


-0.07 


176155-091004A 


8.91 


0.799 


-13.8 


6198 


2.13 


-0.02 


5.7 


9.8 


-0.35 


-0.09 


+0.55 


-0.38 


-0.06 


-0.06 


+0.25 


-0.07 


176155-091004B 


8.96 


0.810 


-13.9 


6256 


2.16 


+0.01 


5.5 


9.8 


-0.36 


-0.09 


+0.54 


-0.39 


-0.10 


-0.06 


+0.27 


-0.07 


176155-091004C 


9.07 


0.834 


-16.5 


6294 


2.20 


+0.00 


5.6 


9.7 


-0.36 


-0.09 


+0.55 


-0.39 


-0.06 


-0.06 


+0.27 


-0.07 


176155-091005A 


9.93 


0.026 


-25.6 


6612 


2.54 


+0.06 


5.7 


9.4 


-0.33 


-0.07 


+0.56 


-0.40 


-0.04 


-0.05 


+0.30 


-0.06 


77 Aql (HD 187929) 












+ 0.09(±0.12) 




-0.20(±0.08) 


+ 0.31(±0.06) 


-0.10(±0.14) 


+ 0.16(±0.07) 


187929-091002A 


6.95 


0.553 


-11.4 


5591 


2.05 


+0.14 


4.3 


10.2 


-0.16 


-0.07 


+0.31 


-0.22 


+0.07 


-0.04 


+0.17 


-0.07 


187929-091003A 


7.94 


0.691 


+3.0 


5735 


2.55 


+0.26 


4.8 


11.7 


-0.07 


-0.05 


+0.24 


-0.18 


+0.04 


-0.03 


+0.29 


-0.07 


187929-091003B 


8.00 


0.700 


+4.3 


5542 


2.01 


+0.15 


4.5 


12.0 


-0.14 


-0.07 


+0.23 


-0.20 


-0.05 


-0.04 


+0.19 


-0.07 


187929-091004A 


8.92 


0.828 


+8.9 


5649 


1.67 


-0.03 


5.1 


16.5 


-0.27 


-0.09 


+0.29 


-0.26 


-0.28 


-0.05 


+0.10 


-0.07 


187929-091004B 


8.97 


0.835 


+ 7.8 


5799 


2.22 


+0.07 


5.7 


16.9 


-0.19 


-0.07 


+0.33 


-0.23 


-0.15 


-0.04 


+0.17 


-0.07 


187929-091004C 


9.08 


0.850 


+5.1 


5654 


1.45 


-0.14 


5.0 


16.6 


-0.32 


-0.11 


+0.31 


-0.29 


-0.32 


-0.05 


+0.04 


-0.07 


187929-091005A 


9.94 


0.970 


-29.4 


6502 


2.61 


+0.14 


5.4 


12.5 


-0.24 


-0.07 


+0.44 


-0.34 


+0.02 


-0.05 


+0.18 


-0.07 
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Table 2. (Continued.) 



Code 


JD 


4> 


^rad 


T eff 


logs 


[Fe] 






[C] 




[N] 


A N 


[O] 


Ao 


[Na] 


A Na 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


S Sge (HD 188727) 












±0 .02(±0 .03) 




-0.25(±0.03) 


±0.39(±0.10) 


-0.10(±0.16) 


±0.15(±0.02) 


188727-091002A 


6.95 


0.702 


+22.3 


5410 


1.59 


-0.01 


4.2 


11.6 


-0.26 


-0.08 


+0.21 


-0.20 


-0.24 


-0.04 


+0.14 


-0.07 


188727-091003A 


7.94 


0.820 


+20.0 


5750 


1.83 


+0.04 


5.2 


14.0 


-0.24 


-0.09 


+0.28 


-0.26 


-0.27 


-0.05 


+0.12 


-0.07 


188727-091003B 


8.01 


0.827 


+ 18.4 


5690 


1.66 


-0.06 


5.5 


14.0 


-0.24 


-0.10 


+0.33 


-0.27 


-0.22 


-0.05 


+0.10 


-0.07 


188727-091003C 


8.10 


0.839 


+ 15.3 


5878 


1.93 


+0.04 


5.6 


13.7 


-0.31 


-0.09 


+0.28 


-0.28 


-0.34 


-0.05 


+0.17 


-0.07 


188727-091004A 


8.94 


0.939 


-9.3 


6176 


2.12 


+0.01 


5.4 


10.8 


-0.26 


-0.09 


+0.44 


-0.36 


+0.02 


-0.07 


+0.15 


-0.07 


188727-091004B 


8.97 


0.943 


-10.4 


6175 


2.15 


+0.04 


5.1 


10.8 


-0.23 


-0.09 


+0.49 


-0.36 


+0.03 


-0.07 


+0.15 


-0.07 


188727-091004C 


9.08 


0.956 


-12.1 


6238 


2.19 


+0.04 


5.0 


10.7 


-0.24 


-0.09 


+0.49 


-0.37 


+0.04 


-0.07 


+0.16 


-0.07 


188727-091005A 


9.94 


0.059 


-12.9 


6109 


1.98 


+0.02 


5.0 


9.7 


-0.21 


-0.10 


+0.49 


-0.38 


+0.07 


-0.08 


+0.14 


-0.08 


188727-091005B 


10.07 


0.074 


-12.2 


6136 


2.08 


+0.06 


4.8 


9.5 


-0.21 


-0.10 


+0.48 


-0.37 


+0.05 


-0.07 


+0.17 


-0.07 



X Cyg (HD 197572) 












±0.09(±0.07) 




-0.30(±0.04) 


±0.49(±0.05) 


±0.01(±0.06) 


+ 0. 21 (±0.04) 


197572-091002A 


7.12 


0.200 


-10.0 


5537 


1.67 


+0.14 


4.2 


9.5 


-0.34 


-0.08 


+0.59 


-0.29 


+0.08 


-0.06 


+0.21 


-0.07 


197572-091003a 


7.95 


0.250 


-4.2 


5349 


1.44 


+0.09 


4.0 


8.3 


-0.33 


-0.09 


+0.53 


-0.24 


+0.11 


-0.05 


+0.19 


-0.07 


197572-091003b 


8.02 


0.254 


-4.3 


5331 


1.32 


+0.06 


3.9 


8.2 


-0.34 


-0.09 


+0.52 


-0.25 


+0.05 


-0.05 


+0.19 


-0.07 


197572-091003c 


8.11 


0.260 


-3.5 


5509 


1.72 


+0.18 


4.0 


8.0 


-0.34 


-0.08 


+0.44 


-0.25 


-0.07 


-0.05 


+0.28 


-0.07 


197572-091004A 


8.99 


0.314 


+0.6 


5059 


0.81 


-0.07 


3.7 


7.3 


-0.32 


-0.10 


+0.47 


-0.22 


+0.01 


-0.05 


+0.13 


-0.05 


197572-091004B 


9.09 


0.320 


+ 1.4 


5167 


1.19 


+0.05 


3.8 


7.2 


-0.28 


-0.09 


+0.50 


-0.21 


+0.01 


-0.04 


+0.17 


-0.06 


197572-091004C 


9.17 


0.325 


+ 2.2 


5243 


1.33 


+0.09 


3.8 


7.0 


-0.29 


-0.09 


+0.48 


-0.22 


-0.03 


-0.04 


+0.20 


-0.07 


197572-091005A 


10.00 


0.375 


+6.7 


5265 


1.54 


+0.14 


3.8 


6.8 


-0.24 


-0.08 


+0.44 


-0.20 


-0.06 


-0.04 


+0.26 


-0.07 


197572-091005B 


10.09 


0.381 


+ 7.0 


5254 


1.56 


+0.14 


3.9 


6.8 


-0.22 


-0.08 


+0.41 


-0.19 


+0.03 


-0.04 


+0.26 


-0.07 


T Vul (HD 198726) 












-0.03(±0.03) 




-0.29(±0.05) 


±0.35(±0.07) 


-0.08 


(±0.07) 


±0.09(±0.02) 


198726-091002A 


6.96 


0.530 


+9.9 


5821 


2.46 


+0.04 


4.5 


10.0 


-0.22 


-0.06 


+0.34 


-0.22 


-0.06 


-0.03 


+0.10 


-0.07 


198726-091003A 


7.96 


0.754 


+ 16.8 


5854 


2.23 


-0.03 


5.2 


12.5 


-0.29 


-0.07 


+0.24 


-0.24 


-0.15 


-0.04 


+0.12 


-0.07 


198726-091003B 


8.03 


0.770 


+ 16.6 


5847 


2.25 


-0.06 


5.4 


12.6 


-0.24 


-0.07 


+0.28 


-0.24 


-0.10 


-0.04 


+0.09 


-0.07 


198726-091003C 


8.12 


0.792 


+ 15.0 


5917 


2.20 


-0.03 


4.9 


12.7 


-0.33 


-0.08 


+0.27 


-0.27 


-0.23 


-0.04 


+0.09 


-0.07 


198726-091004A 


8.95 


0.979 


-17.5 


6558 


2.39 


-0.09 


4.7 


11.6 


-0.35 


-0.08 


+0.37 


-0.39 


-0.04 


-0.06 


+0.07 


-0.07 


198726-091004B 


9.00 


0.990 


-17.6 


6605 


2.48 


-0.05 


4.4 


11.5 


-0.33 


-0.08 


+0.38 


-0.40 


-0.07 


-0.06 


+0.08 


-0.07 


198726-091004C 


9.10 


0.013 


-18.3 


6562 


2.48 


-0.04 


4.6 


11.2 


-0.33 


-0.08 


+0.39 


-0.39 


-0.04 


-0.06 


+0.10 


-0.07 


198726-091005A 


9.98 


0.211 


-7.5 


6036 


2.38 


-0.02 


4.3 


8.9 


-0.22 


-0.07 


+0.46 


-0.30 


+0.02 


-0.05 


+0.07 


-0.08 


198726-091005B 


10.10 


0.237 


-6.8 


6010 


2.31 


-0.02 


4.3 


8.7 


-0.28 


-0.08 


+0.39 


-0.30 


-0.02 


-0.05 


+0.06 


-0.08 



DT Cyg (HD 201078) 










±0.08(±0.06) 




-0.19 


(±0.05) 


±0.49(±0.03) 


±0.05(±0.02) 


±0.25(±0.03) 


201078-091003A 


7.97 


0.820 


+ 2.6 


6245 


2.43 


-0.02 


4.1 


10.1 


-0.24 


-0.08 


+0.44 


-0.34 


+0.04 


-0.06 


+0.18 


-0.07 


201078-091003B 


8.04 


0.850 


+ 1.0 


6341 


2.58 


-0.01 


4.3 


10.4 


-0.24 


-0.07 


+0.44 


-0.34 


+0.04 


-0.05 


+0.23 


-0.07 


201078-091003C 


8.13 


0.886 


-0.5 


6448 


2.82 


+0.10 


4.2 


10.7 


-0.17 


-0.07 


+0.49 


-0.33 


+0.03 


-0.05 


+0.26 


-0.07 


201078-091003D 


8.19 


0.911 


-2.0 


6502 


2.92 


+0.14 


4.5 


11.1 


-0.08 


-0.06 


+0.47 


-0.32 


+0.07 


-0.05 


+0.28 


-0.07 


201078-091004A 


9.01 


0.239 


-4.7 


6386 


2.75 


+0.11 


3.9 


8.3 


-0.19 


-0.07 


+0.51 


-0.34 


+0.07 


-0.05 


+0.27 


-0.07 


201078-091004B 


9.11 


0.278 


-3.5 


6292 


2.65 


+0.08 


3.8 


8.1 


-0.19 


-0.07 


+0.52 


-0.34 


+0.08 


-0.05 


+0.25 


-0.07 


201078-091004C 


9.18 


0.307 


-3.0 


6421 


3.04 


+0.18 


4.0 


7.8 


-0.15 


-0.06 


+0.55 


-0.30 


+0.08 


-0.04 


+0.30 


-0.07 


201078-091005A 


10.02 


0.642 


+ 5.2 


6181 


2.71 


+0.07 


4.1 


8.6 


-0.20 


-0.06 


+0.49 


-0.29 


+0.03 


-0.04 


+0.24 


-0.07 


201078-091005B 


10.11 


0.680 


+ 5.0 


6155 


2.52 


+0.04 


4.0 


9.1 


-0.26 


-0.07 


+0.51 


-0.32 


+0.02 


-0.05 


+0.21 


-0.07 


V1334 Cyg (HD 203156) 










±0.01(±0.06) 




-0.24 


(±0.03) 


±0.44(±0.03) 


±0.03(±0.03) 


±0.11(±0.04) 


203156-091003A 


8.06 


0.886 


-22.4 


6286 


2.20 


-0.06 


4.0 


14.7 


-0.26 


-0.10 


+0.48 


-0.40 


+0.05 


-0.07 


+0.03 


-0.07 


203156-091003B 


8.14 


0.911 


-22.6 


6496 


2.60 


+0.03 


4.3 


14.7 


-0.26 


-0.08 


+0.46 


-0.38 


+0.04 


-0.06 


+0.13 


-0.07 


203156-091003C 


8.20 


0.927 


-22.6 


6503 


2.65 


+0.05 


4.4 


14.6 


-0.28 


-0.07 


+0.45 


-0.37 


+0.07 


-0.06 


+0.12 


-0.07 


203156-091004A 


8.93 


0.145 


-19.0 


6260 


2.33 


-0.03 


4.0 


14.4 


-0.23 


-0.09 


+0.43 


-0.36 


+0.04 


-0.06 


+0.09 


-0.07 


203156-091004B 


9.03 


0.176 


-19.2 


6231 


2.21 


-0.01 


3.7 


14.5 


-0.25 


-0.10 


+0.42 


-0.38 


+0.00 


-0.07 


+0.09 


-0.08 


203156-091004C 


9.12 


0.204 


-17.9 


6388 


2.66 


+0.07 


4.1 


14.6 


-0.20 


-0.07 


+0.45 


-0.34 


+0.02 


-0.05 


+0.15 


-0.07 


203156-091004D 


9.19 


0.226 


-17.1 


6246 


2.30 


-0.03 


4.0 


14.4 


-0.27 


-0.09 


+0.43 


-0.37 


+0.01 


-0.06 


+0.10 


-0.08 


203156-091005A 


10.04 


0.479 


-12.2 


6249 


2.26 


-0.04 


3.9 


14.8 


-0.22 


-0.10 


+0.39 


-0.37 


-0.03 


-0.06 


+0.08 


-0.08 


203156-091005C 


10.18 


0.523 


-11.5 


6499 


2.91 


+0.11 


4.5 


15.0 


-0.18 


-0.06 


+0.47 


-0.33 


+0.04 


-0.05 


+0.18 


-0.07 



S Cep (HD 213306) 




















±0.04(±0.05) 






-0 


.22 


(±0. 


05) 


±0.39(±0. 


07) 


±0 


.02(±0.06) 


±0.17(±0. 


04) 


213306-091002A 


6. 


.98 





.502 


-14. 


,1 


5706 


2 


.19 


+o. 


.07 


3.9 


7. 


.0 


-0. 


23 


-0 


.07 


+0 


.40 


-0, 


,24 


-0, 


.01 


-0.04 


+0.15 


-0, 


,07 


213306-091002B 


7. 


.22 





.545 


-11. 


.5 


5713 


2. 


.31 


+o. 


.09 


4.2 


7, 


.6 


-0. 


18 


-0 


.06 


+0 


.41 


-0. 


,22 


+o. 


.04 


-0.04 


+0.17 


-0. 


,07 


213306-091003A 


7. 


.98 





.687 


-4. 


.1 


5682 


2, 


.06 


+o. 


.11 


4.5 


9 


.8 


-0. 


24 


-0 


.07 


+o 


.29 


-0. 


.23 


-0, 


.03 


-0.04 


+0.21 


-0. 


.07 


213306-091003B 


8. 


.07 





.705 


-3 


.8 


5661 


2, 


.12 


+o. 


.10 


4.6 


10, 


.2 


-0. 


20 


-0 


.07 


+o 


.33 


-0. 


,22 


+o. 


.02 


-0.04 


+0.19 


-0. 


,07 


213306-091003C 


8. 


.15 





.720 


-1 


.9 


5666 


2 


.14 


+o. 


.09 


4.8 


10, 


.5 


-0. 


17 


-0 


.07 


+o 


.34 


-0. 


,22 


+o. 


.03 


-0.04 


+0.18 


-0. 


.07 


213306-091003D 


8. 


.21 





.730 


-1. 


,5 


5625 


1. 


.98 


+o. 


,01 


4.8 


10. 


.8 


-0. 


18 


-0 


.07 


+o 


.33 


-0, 


,23 


+o. 


.04 


-0.04 


+0.17 


-0, 


.07 


213306-091003E 


8. 


.26 





.739 


-1 


.6 


5633 


1. 


.92 


+o. 


.03 


4.8 


11. 


.0 


-0. 


17 


-0 


.08 


+o 


.30 


-0, 


,23 


+o. 


.01 


-0.04 


+0.17 


-0. 


,07 


213306-091003F 


8. 


.30 





.747 


-0 


.1 


5658 


1. 


.99 


+o. 


.03 


4.8 


11 


.2 


-0. 


14 


-0 


.08 


+o 


.32 


-0. 


.23 


+o. 


.01 


-0.04 


+0.18 


-0. 


,07 


213306-091004A 


8. 


.96 





.871 


-11, 


.2 


6330 


2. 


.66 


+o. 


.10 


6.0 


11 


.7 


-0. 


28 


-0 


.06 


+o 


.32 


-0, 


,28 


-0 


.09 


-0.04 


+0.24 


-0. 


,07 


213306-091004B 


9. 


.05 





.886 


-17. 


,7 


6319 


2. 


.52 


+o. 


.03 


5.8 


10. 


.9 


-0. 


,27 


-0 


.07 


+o 


.40 


-0. 


,31 


-0, 


.01 


-0.05 


+0.20 


-0. 


.07 


213306-091004C 


9. 


.13 





.902 


-22 


,7 


6323 


2. 


.50 


+o. 


.01 


5.6 


10 


.4 


-0. 


22 


-0 


.07 


+o 


.48 


-0. 


,33 


+o. 


.07 


-0.05 


+0.16 


-0, 


.07 


213306-091004D 


9. 


.21 





.916 


-27. 


,1 


6545 


2. 


.65 


+o. 


,05 


5.6 


10. 


.3 


-0. 


26 


-0 


.07 


+o 


.43 


-0. 


,34 


+o. 


.02 


-0.05 


+0.20 


-0. 


.07 


213306-091004E 


9. 


.27 





.929 


-31, 


.1 


6663 


2. 


.51 


-0 


.03 


5.3 


10. 


.4 


-0. 


33 


-0 


.07 


+o 


.38 


-0. 


.38 


-0. 


.07 


-0.05 


+0.18 


-0. 


.06 


213306-091005A 


9. 


.96 





.058 


-36 


.9 


6515 


2. 


.31 


-0 


.09 


5.0 


9. 


.3 


-0. 


,27 


-0 


.09 


+o 


.44 


-0. 


,41 


-0 


.02 


-0.06 


+0.08 


-0. 


,07 


213306-091005B 


10. 


.21 





.104 


-34, 


.7 


6427 


2. 


.44 


+o. 


.08 


4.5 


9. 


.4 


-0. 


22 


-0 


.08 


+o 


.47 


-0. 


,38 


+o. 


.07 


-0.06 


+0.19 


-0. 


.07 


213306-091005C 


10. 


.27 





.114 


-34. 


,3 


6189 


2. 


.12 


+o. 


,05 


4.4 


9. 


.2 


-0. 


21 


-0 


.10 


+o 


.56 


-0, 


,40 


+o. 


.18 


-0.08 


+0.10 


-0. 


.08 



Column (1) — spectrum code (??????-yymmdd#), where ?????? is the HD number, yymmdd is the observation date (UT), and # denotes the spectrum turn of 
a star for the day (A ■ ■ ■ 1st, B ■ ■ ■ 2nd, etc.). Column (2) — observational time in the Julian day expressed as JD — 24455100. Column (3) — pulsation phase. 
Column (4) — heliocentric radial velocity (in km s — 1 ). In Columns (5)— (9) are given the results of atmospheric parameters, T e ff (effective temperature, in K), 
log g (surface gravity, in cm s — 2 ), [Fe/H] (logarithmic Fe abundance relative to the Sun), £ (microturbulence, in km s — 1 ), and vy[ (macrobroadening velocity 
in km s — 1 derived from 6143—6163 A fitting), respectively. The last 8 Columns present the final results of [X/H] (logarithmic abundance of element X relative 
to the Sun) and A x (non-LTE correction) for C (10, 11), N (12, 13), O (14, 15), and Na (16, 17). The mean abundance (([X/H])) averaged over each of the 
phases along with the standard deviation (cr x ) are also given at the first line of each section (expressed in italic). 



C, N, 0, and Na Abundances of Cepheid Variables 



Table 3. Adopted atomic data of important spectral lines. 



Species 


A (A) 


Xlow (eV) 




Gammar 


Gammas 


Gammaw 


Remark 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 








[6143-6163 A fitting] 






Si i 


61 4^ 01 6 


5.616 


—0.820 


(7.77) 


(—4.45) 


( — 7 0^) 




Fe II 


fil 47 741 

U14 1 . Ill 


3.889 


—2.721 


8.53 


—6.53 


—7.88 




Fe II 


(S14Q 
\j ±'-tty . 


3.889 


—2.724 


8.53 


—6.53 


—7.88 




Fe II 


fil ^0 0Q& 


3.221 


—4.754 


8.54 


—6.54 


—7.91 




Fe I 


61 ^1 61 7 

UJ.tJ-L.U-L 1 


2.176 


—3.299 


8.19 


—6.20 


—7.82 




Na I 


fil ^4 99fi 


2.102 


— 1.560 


7.85 


—4.39 


(—7.29) 


Npi t fil ^4 

1 N ex 1 U-LO^t 


Si i 


61 ^ 1 S4 


5.619 


-0.400 


(7.77) 


(-4.45) 


(—7.05) 




Si i 


fil ^ fiQ^ 


5.619 


-1.690 


(7.77) 


(-4.45) 


(—7.05) 




O i 


61 Q61 


10.740 


-1.401 


7.60 


-3.96 


(—7.23) 


O t 61 "^-8 


O i 


fil ^ Q71 

U-L tJO. c/ 1 J. 


10.740 


-1.051 


7.61 


-3.96 


(—7.23) 


O t fil 


O i 


6155.989 


10.740 


-1.161 


7.61 


-3.96 


(-7.23) 


O I 6155-8 


Ca i 


6156.023 


2.521 


-2.200 


7.49 


-4.69 


-7.50 




O i 


6156.737 


10.740 


-1.521 


7.61 


-3.96 


(-7.23) 


O I 6155-8 


O i 


6156.755 


10.740 


-0.931 


7.61 


-3.96 


(-7.23) 


O I 6155-8 


O i 


6156.778 


10.740 


-0.731 


7.62 


-3.96 


(-7.23) 


O I 6155-8 


Fe i 


6157.725 


4.076 


-1.260 


7.70 


-6.06 


-7.84 




O i 


6158.149 


10.741 


-1.891 


7.62 


-3.96 


(-7.23) 


O I 6155-8 


O i 


6158.172 


10.741 


-1.031 


7.62 


-3.96 


(-7.23) 


O I 6155-8 


O i 


6158.187 


10.741 


-0.441 


7.61 


-3.96 


(-7.23) 


O I 6155-8 


Fe I 


6159.368 


4.607 


-1.970 


8.28 


-4.61 


-7.77 




Na I 


6160.747 


2.104 


-1.260 


7.85 


-4.39 


(-7.29) 


Na I 6161 


Ca I 


6161 297 


2.523 


-1.020 


7.49 


-4.69 


—7.50 




Ca I 


6162.173 


1.899 


+0.100 


7.82 


-5.07 


-7.59 










[7110-7121 A fitting] 






Ni I 


7110.892 


1.935 


-2.880 


7.72 


-6.30 


-7.85 


(adjusted gf) 


C i 


7111.472 


8.640 


-1.240 


(7.64) 


(-5.02) 


(-7.24) 


C I 7111, (adjusted gf) 


C i 


7113.178 


8.647 


-0.800 


(7.64) 


(-5.02) 


(-7.23) 


C I 7113, (adjusted gf) 


C i 


7115.172 


8.643 


-0.960 


(7.64) 


(-5.02) 


(-7.23) 


C I 7115, (adjusted gf) 


C i 


7115.182 


8.640 


-1.550 


(7.64) 


(-5.02) 


(-7.24) 


C I 7115 


C i 


7116.991 


8.647 


-0.910 


(7.64) 


(-5.02) 


(-7.23) 


C I 7116 


Fe i 


7118.119 


5.009 


-1.390 


8.69 


-5.22 


-7.75 


(adjusted gf) 


C i 


7119.656 


8.643 


-1.130 


(7.64) 


(-5.02) 


(-7.23) 


C I 7119, (adjusted gf) 


Fe I 


7120.022 


4.558 


-1.911 


8.55 


-5.28 


-7.68 


(adjusted gf) 








[8677-8697 A 


. fitting] 






S i 


Rf\7R Q97 

OU 1 O. c/Zi 1 


7.867 


-1.000 


7.61 


-4.41 






S i 


8679.620 


7.867 


-0.410 


7.61 


-4.41 


(-7.30) 




Si i 


8680.080 


5.862 


-1.000 


(7.47) 


(-5.00) 


(-7.23) 




N i 


8680.282 


10.336 


+0.236 


8.62 


-5.51 


(-7.63) 


N I 8680 


S i 


8680.411 


7.867 


-0.210 


7.61 


-4.41 


(-7.30) 




N i 


8683.403 


10.330 


-0.045 


8.62 


-5.51 


(-7.64) 


N I 8683 


N I 


8686.149 


10.326 


-0.448 


8.62 


-5.51 


(-7.64) 


N I 8686 


Si i 


8686.352 


6.206 


-0.700 


(7.47) 


(-5.00) 


(-7.23) 


(adjusted gf) 


Fe I 


8688.621 


2.176 


-1.212 


7.28 


-6.24 


-7.85 




Fe I 


8689.857 


5.105 


-1.948 


8.07 


-5.77 


-7.73 


(adjusted gf) 


Si I 


8690.061 


5.613 


-1.980 


(7.47) 


(-5.28) 


(-7.32) 




S I 


8693.137 


7.870 


-1.370 


7.62 


-4.41 


(-7.30) 




S i 


8693.931 


7.870 


-0.510 


7.62 


-4.41 


(-7.30) 




S i 


8694.626 


7.870 


+0.080 


7.62 


-4.41 


(-7.30) 





Note. All data are were taken from Kurucz & Bell's (1995) compilation as far as available, though empirically adjusted 
values (solar gf values) were used for several lines remarked as "adjusted" in Column (8). The meanings of Columns (l)-(4) 
are self-explanatory. In Columns (5)-(7) are given the damping parameters in the c.g.s. unit: 

Gammar is the radiation damping constant (s _1 ), log7 ra d. Gammas is the Stark damping width per electron density (cm -3 ) 
at 10 4 K, log(7e/Af e ). Gammaw is the van der Waals damping width per hydrogen density (cm -3 ) at 10 4 K, log(7 w /A^H)- 
Note that the values in parentheses are the default damping parameters computed within the Kurucz's (1993) WIDTH9 
program (cf. Leusin & Topil'skaya 1987), because the damping data for these lines were unavailable in Kurucz & Bell's 
(1995) database. 
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Figure 1. Line-by-line Fe abundances (corresponding to the finally established atmospheric 
parameters), plotted against equivalent widths (left panels) and lower excitation potentials 
(right panels), for the case of 5 Cep (213306-091002A; the first spectrum for this star). The 
upper two panels (a, b) are for Case (A) where lines of w < 200 mA (w is the reduced 
equivalent width: w = W\ • (5000/A)) were equally used for both Fe I and Fe II lines, while 
the lower two panels (c, d) correspond to Case (B) where lines were screened with the criterion 
of w < 30 mA (Fe i) and w < 200 mA (Fe n) (i.e., only weak Fe I lines were used). The 
results for Fe I and Fe II lines are distinguished by filled and open symbols, respectively, and 
the average Fe abundance is indicated by the horizontal dashed line. 
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Figure 2. Fe abundance vs. equivalent width relations corresponding to the finally established atmospheric pa- 
rameters of T e ff, logg, and vt for each of the 122 spectra, being arranged according to the time sequence in the 
downward direction for each star (just as in Table 2). The filled and open symbols correspond to Fe I and Fe II 
lines, respectively. The results for each stars are shown relative to the mean abundances indicated by the horizontal 
dotted lines, and vertically shifted by 1.0 relative to the adjacent ones. 
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Figure 3. Fe abundance vs. lower excitation potential relation corresponding to the finally established atmosptu 
parameters for each of the 122 spectra. Otherwise, the same as in Fig. 2. 
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Figure 4. Synthetic spectrum fitting in the 6143-6163 A region accomplished by varying the abundances of O, Na, 
Si, Ca, and Fe, along with the macrobroadening velocity (vm) an d the wavelength shift (radial velocity). The best-fit 
theoretical spectra are shown by solid lines, while the observed data are plotted by symbols, where the wavelength 
scale of the stellar spectrum has been adjusted to the laboratory frame. In each panel, the spectra are arranged 
according to the time sequence in the downward direction for each star (just as in Table 2), and each spectrum is 
vertically shifted by 0.25 relative to the adjacent one. The important lines (cf. Table 3) are identified in the lowest 
spectrum of the left panel. 




Figure 5. Synthetic spectrum fitting in the 7110-7121 A region accomplished by varying the abundances of C, Fe, 
and Ni. Each spectrum is vertically shifted by 0.15 relative to the adjacent one. Otherwise, the same as in Fig. 4. 
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Figure 6. Synthetic spectrum fitting in the 8677-8697 A region accomplished by varying the abundances of N, Si, 
S, and Fe. Each spectrum is vertically shifted by 0.25 relative to the adjacent one. Otherwise, the same as in Fig. 4. 
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Figure 7. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 17 spectra of SU Cas, plotted 
against the pulsation phase. Shown in these 7 panels are (from top to bot- 
tom): (a) heliocentric radial velocity (km s _1 ), (b) effective temperature (K), 
(c) logarithmic surface gravity (cm s -2 ), (d) microturbulence (km s _1 ), (e) 
macrobroadening velocity (km s _1 ), (f) [Fe/H] (dex), and (g) [X/H] (dex) 
where X is C (filled triangles), N (filled circles), O (open circles), and Na 
(greek crosses). 
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Figure 8. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 12 spectra of SZ Tau, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 9. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 10 spectra of RT Aur, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 10. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 8 spectra of £ Gem, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 11. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 7 spectra of FF Aql, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 12. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 7 spectra of 77 Aql, plotted against 
the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 13. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 9 spectra of S Sge, plotted against 
the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 14. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 9 spectra of X Cyg, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 15. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 9 spectra of T Vul, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 16. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 9 spectra of DT Cyg, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 17. Results of the radial velocity, atmospheric parameters, and ele- 
mental abundances derived from each of the 9 spectra of V1334 Cyg, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 18. Results of the radial velocity, atmospheric parameters, and el- 
emental abundances derived from each of the 16 spectra of 5 Cep, plotted 
against the pulsation phase. Otherwise, the same as in Fig. 7. 
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Figure 19. Phase (</>)-dependence of the atmospheric parameters of 5 Cep (HD 213306), 
where our results and the previous results of other similar "spectroscopic" determinations 
are overplotted with different symbols: Filled circles — this study, open circles — Kovtyukh 
& Andrievsky (1999) (spectroscopic T e ^ and non-standard results for logg, £,and [Fe/H]), 
open triangles — Fry & Carney (1997), double circle — Kovtyukh et al. (1994) (note that 
their T e ^ was determined from Ha and colors), Greek crosses — Luck Sz Lambert (1985) 
(spectroscopic logg), diamonds — Andrievsky et al. (2005). (a) T e ^ vs. <f>, (b) logg vs. </>, (c) 
£ vs. </>, and (d) [Fe/H] vs. <j>. 



C, N, 0, and Na Abundances of Cepheid Variables 



7000h 

g 

^ 6000 



5000- 




0.2 0.4 0.6 0.8 1 1.2 
log P (day) 



3 2- 



1 I 1 I 1 I 1 I 



3 _ i a + . 



O 



1- 



(b) 



J i I i I i L 



k2 0.4 0.6 0.8 1 
log P (day) 



1.2 



"i — 1 — i — 1 — i — 1 — i — 1 — r 



° i 



i I i I i I i_ 



_L 



0.2 0.4 0.6 0.8 1 1.2 
log P (day) 



15- 



1 l 1 I 1 I 1 I 



8 1 » . x • 



5 -(d), 







0.2 0.4 0.6 0.8 1 
log P (day) 



1.2 



0.5 



^3 



-0.5- 



1 l 1 l 1 l 1 l 1 l 



H M'i! !|i U 



(e). 



i ■ i 



0.2 0.4 0.6 0.8 1 1.2 
log P (day) 



Figure 20. Atmospheric parameters derived from each of the 122 spectra, plotted against 
the pulsation period (logP): (a) effective temperature, (b) logarithmic surface gravity, (c) 
microturbulence, (d) macrobroadening velocity, and (e) [Fe/H]. In panel (b), an approximate 
relation for dynamical logg given by equation (8) is also depicted by a dashed line. The 
results for each star are distinguished by the symbol type: SU Cas — open circles, DT Cyg 

— filled triangles, SZ Tau — open triangles, V1334 Cyg — Greek crosses (+), RT Aur — 
filled squares, T Vul — open downward triangles, FF Aql — St. Andrew's crosses (x), S Cep 

— filled circles, 77 Aql — open diamonds, S Sge — filled downward triangles, £ Gem — open 
squares, X Cyg — filled diamonds. 
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Figure 21. ([X/H]) (mean logarithmic abundance of element X relative to the Sun averaged 
over different pulsation phases) for each star, plotted against the pulsation period, based on 
the data (expressed in italic) given in the first line of each section in Table 2: (a) ([C/H]), 
(b) ([N/H]>, (c) ([O/H]}, (d) ([Na/H]), and (e) ([Fe/H]). The same meanings of the symbols 
as in Fig. 20, while the attached error bar indicate the extent of the standard deviation (<x). 
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Figure 22. Panels (a)-(d) display the correlations between the C, N, O, and Na abun- 
dances (averaged over different phases): (a) ([N/H]> vs. ([C/H]>, (b) ([Na/H]) vs. ([O/H]), 
(c) (loge(N)) vs. (loge(C+0)>, and (d) (loge(N)) vs. (loge(C)) [(loge(X)) = ([X/H]> + 
loge©]. In panels (c) and (d), the expected relation when e(C)+e(N)+e(0) is conserved at 
the solar value is also shown. The solid and dashed line in panel (d) correspond to the cases 
of [O/H] = 0.0 and [O/H] = -0.1, respectively. Further, the N-to-C abundance ratio ([N/C] 
= [N/H]-[C/H]) and the sum of C+N+O abundances ((log e(C+N+0)>) for each star are 
plotted against logP in panels (e) and (f), respectively. Otherwise, the same as in Fig. 20. 



